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Abstract. 
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on spectrum generating algebras. 
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1. Introduction: Darboux intertwining, Schrodinger factorization, Witten 
SUSY mechanics in one basket 

The concept of Supersymmetric Quantum Mechanics (SUSY QM) initially associated 
to the + 1 dimensional SUSY Field Theory [l] aimed at a simplified analysis of difficult 
problems of mult i- dimensional QFT such as spontaneous SUSY breaking, vacuum 
properties beyond perturbation theory etc. In addition, the realization of different 

SUSY algebras in Quantum Physics turned out to be easily achieved in certain QM 
models [U |3] . 

Soon after its formulation SUSY QM was well identiffed with the Quantum 
Mechanics of isospectral systems described by Hamiltonians with almost coincident 
energy spectra 111^^0] • SUSY manifested itself through intertwining Darboux [TT] 
transformations between isospectral partners, used before in Mathematics 52]^PE]- The 
latter property presumably gave E. Schrodinger ^Tj a hint to factorize the corresponding 
Hamiltonians into a product of simplest, ffrst-order Darboux operators ^H]- Altogether 
the differential realization of SUSY in QM stimulated its application for the spectral 
design, i.e. for preparation of quantum potential systems with pre-planned energy 



spectra [IHl^jSni and scattering data [lOl^lSn] or potential profiles jHIl E21 in a 



constructive way. For the spectral design-kit the non-linear supersymmetry in its 
differential realization |SS]^|nni has become one of the most efficient tools to build 
various isospectral quantum systems with desired features. By now, several books 
and reviews devoted to certain achievements and diverse applications of 

SUSY QM approach illuminated at length many of the above mentioned trends. In 
contradistinction, the role and the structure of nonlinear SUSY QM just developed in 
the last decade and has not been yet surveyed. 

The interplay between the algebraic and differential properties of nonlinear 
supersymmetry, in the spectral design, is guiding our present work: we are going to 
clarify the benefits of algebraic SUSY approach to the old Darboux-Crum method, to 
develop, with these tools, the complete classification of differential realizations of non- 
linear SUSY algebras and to give insight to intrinsic links between SUSY isospectrality 
and hidden symmetries in particular quantum systems. 

In what follows we restrict ourselves to the full analysis in the one-dimensional, one- 
component, stationary QM elaborating few important examples for the non-stationary 
Schrodinger equation and leave complex, matrix, multidimensional and/or relativistic 
(Dirac, Klein-Gordon etc.) equations without any detailed comments. 

Let us start the retrospection of SUSY QM: consider two one-dimensional 
Schrodinger Hamiltonians defined on the line and assemble them into a matrix 
Super-Hamiltonian, 



with non-singular real potentials. These Hamiltonians h"^ are supposed to have (almost) 
the same energy levels for bound states and/or the same spectral densities for continuum 
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parts of spectra. Furthermore, assume that their isospectral connection is provided by 
intertwining with the Crum-Darboux fn\ IT^ transformation operators q^, 

h^QN = Qn^-, q^h^ = h~qj^, (2) 
where iV-th order differential operators 

N 

q% = Y.^t{x)d\ = const ^ {Tl)"" . (3) 

fc=0 

The conventional, linear M = 1 SUSY QM in the fermion number representation [2j is 
implemented by nilpotent supercharges Qi, Q[ of first order in derivatives built from a 
real super-potential 

qf^^d + x{xY =^ ^^=(o f )' ^?=(^0' = O' 

where ^ stands for the operation of hermitian conjugation (as well as we employ the 
differential operators with real coefficients the hermitian conjugation is equivalent to 
the operation of transposition, = Q*) . 

The intertwining relations introduced in result in the Supersymmetry for a 
Super-Hamiltonian H, 

[H,Q,] = [H,Q{] = 0. (5) 

The SUSY algebra is completed by the appropriate decomposition of the Super- 
Hamiltonian, 

H = {Q,, Q\} ^ h+ = q+q- = + - x'; = q^qt = -d' + x' + x' , (6) 

which is in line with the Schrodinger one-step factorization ^3 ^j- The notation 
x' = dx/dx has been employed. 

At this stage, the super-potential x is generated by zero-energy solutions of the 
Schrodinger equations (equivalently, zero-modes of supercharges Qi,Qi), 

h^'PTo) = = Qi'P'ior <^ro)(^) = i'Pto)^^)) ^ = exp (^- J dyx{y)j ■ (7) 

If x{x) is a non-singular function then the zero-modes 0^-, are nodeless. This leads to 
a non-negative physical spectrum in agreement with the SUSY algebra ©, 

= iQii^lQii^) + {QUIQU) > o, (8) 

for any L^-normalizable, smooth wave function ijj[x). 

Several options exist for the choice of zero-modes of supercharges Qi,Qi- If one of 
the zero-modes is normalizable then it becomes a ground state wave function of the 
Super-Hamiltonian H {i.e. of or h^). But another one remains non- normalizable 
due to Eq.((Zj). Thus either q^ or q^ deletes the ground state level of /i^ or . When 
keeping in mind the spectral design program one can also interpret it conversely: if q~ 
deletes the lowest level of h'^ converting it into h" then g+ creates a new level for h~ 
transforming it into h'^. 
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Another option is realized by the non-normahzable nodeless functions when 
none of them belongs to the physical spectrum of Hamiltonians h^. In this case the 
entire sets of physical eigenstates of the both Hamiltonians are put into the one-to-one 
correspondence by intertwining relations 

h^^^ = E^%- E>0- ^| = -^g,T^|, (9) 

and such Hamiltonians are strictly isospectral. In the SUSY vocabulary it is the case 
of "spontaneous" SUSY breaking as the lowest ground state of the Super-Hamiltonian 
H is degenerate. 

The previous analysis has been based not only on the intertwining relations (0) but 
also on the factorization However there is no such factorization for higher-order 
intertwining operators Q- What do we have for the latter ones instead? 



2. From the ladder of SUSY's via Parasupersymmetry toward Polynomial 
SUSY 

Let us proceed by recursion and discover different levels of isospectrality: from a 
simple Darboux transformation to a ladder or a dressing chain made of several simple 
Darboux steps. Actually the whole variety of elementary building blocks for spectral 
design can be well developed within the class of transformation operators qf of second- 
order in derivatives. One has to select the operators (jH)) with nonsingular coefficient 
functions which produce a nonsingular potential V2 after intertwining (0) with the 
smooth initial potential Vi. 

First of all, to produce the required transformation operators the two different linear 
SUSY systems may be "glued". Indeed, consider two Super-Hamiltonians Hi, i = 1,2, 
Eq. ((H), respectively two supercharges Qi with super-potentials Xi and supercharge 
components rf = + Xi- Let us identify two elements of Super-Hamiltonians, 

h^ = ht + X; xi + x'i = xl-X2 + >^ (10) 

with e[^^ > a > —E^^ where E^^ and Ef"* are ground state energies for hi and /i^ 
respectively. Evidently the constant shift of the Super-Hamiltonian H2 —>■ H2 + X does 
not break or change the supersymmetry. 

After such a gluing the chain of intertwining relations Q can be assembled into 
the supersymmetry transformation [Hps,Qps] = [Hps,Q'jjg] = of the combined Super- 
Hamiltonian Hps and the joint supercharges Qpg, Qls, 
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where we have shifted both Super-Hamiltonians symmetrically to simplify the further 
algebra. However these supercharges Qps,Q\,s are not anymore nilpotent of order two 
and therefore do not mimic the Pauli principle. Still they are nilpotent of order 
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three, Q^^ = (QIs) — 0' therefore the states carrying these charges obey the 
para-statistical principle. Thus we deal now with the Parasupersymmetry |l()4j - pTin] . 
Furthermore the closure of superalgebra is not anymore given by Eq.®- The lowest- 
order relation between the Super-Hamiltonian and the para-supercharges is trilinear, 

QpsQps ~^ QpsQps ~^ QpsQpsQps — "^HpsQps- (12) 

This quantum system reveals triple degeneracy of levels with the possible exception 
for two lowest states. We draw also the reader's attention to a possibility to treat both 
intermediate Hamiltonians and /i^ + A separately, in spite of their identification, thus 
doubling of the Hilbert spaces spanned by their eigenfunctions. It leads to a model of 
"Weak Supersymmetry" with quadruple level degeneracy jllUI which may be susceptible 
to prolongation onto a weak-SUSY field theory. 

Going back to the spectral design purposes, the para-supersymmetric dynamics 
contains redundant information, namely, about the intermediate Hamiltonian = 
h2 + X. One is, in fact, interested in the final Hamiltonian /ig only as produced from 
the initial one, hf by means of a second-order Darboux transformation, 

/i+r+r+ = r+ /i^r+ = r+(/i+ + A)r+ = r+ r+(/i2 + A); r^r^h^ = {h:^ + Xy^^r^ . (13) 

Let us make a shortcut and define the two isospectral components /i^, 

/i+ = /i+ + Ai = r+rf + Ai; = + X2 = r^^r^ + X2] (14) 
r^r^ + Ai = r^r^ + A2; (15) 

for the generalized Super-Hamiltonian (PJ) where we have employed a more general 
energy reference (shift by arbitrary Ai 2). Evidently, A = A2 — Ai . Then the intertwining 
relations (j21) are identical to Eq. ()13|) with = r^^rj and the supersymmetry 
[H, Q2] = [H, QI] = is generated by the conserved supercharges, 

Q2=( ° i), Ql=(QlY = o. (16) 



^0 ^ 

In place of Eq.®, because of (fT^ the algebraic closure is given by, 

{Q2, Ql} = ( '''''''' _ _\ A={H-X,){H- A2). (17) 
\ r2rir{ri J 

Thus we have obtained the second-order Polynomial SUSY algebra [SHI- [20] as a concise 
form of isospectral deformation of a potential system accomplished by a ladder [Uj-pU], 
^H]^|2S] or a dressing chain j2n]^|S21 of a couple of one-step Darboux transformations 
or, equivalently, by a second-order Crum-Darboux intertwining ^21 EZ] or by 
a blocking of two linear SUSY with partial overlapping of Super-Hamiltonians [53] 
("weak SUSY" M), or by a tower of para-SUSY transformations [TTn]-[TTlK]. 

In the modern SUSY vocabulary there are several synonyms for the higher-order 
SUSY algebra: originally it was named as a Polynomial (or Higher-derivative) one 
[H31IS1], recently the title of A^-fold SUSY has been suggested [73] and, at last, a more 
general term of Nonlinear SUSY has been used [HD] with a certain reference to nonlinear 
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SUSY algebra arising in the conformal QM In what follows we will combine the 

first name and the last one depending on the structure of a superalgebra |85j . 

This Polynomial SUSY keeps track of essential spectral characteristics of the 
second-order SUSY (Crum-Darboux transformations) . Indeed, the zero-modes of 
intertwining operators qf or, equivalently, the zero-modes of the hermitian supercharges 
Qt = Q2 + Qh Q2 = ^(^2 ~ Q2), form the basis of a finite-dimensional representation 
of the Super-Hamiltonian, 

q^<pf{x) = = qih^<Pt{xy, 2=1,2; = (18) 

i=i 

due to intertwining relations ((21), (fT!?|) . In terms of these Hamiltonian projections - 
constant matrices , the SUSY algebra closure takes the polynomial form [H3] (see 
also (73j). 



{Q2,qI} 



det 



EI - 



= det 

E=H 



EI-S- 



^^^^V,{H). (19) 



Thus both matrices have the same set of eigenvalues which for the ladder 
construction (fT7|) consists of Ai,A2 . As the zero-mode set is not uniquely derived 
from ()18|) the matrices are not necessarily diagonal. For instance, the equation 
rj*"r2"0+(x) = has one zero- mode obeying r2'02^(x) = and another one obeying 
'^I 'Pi = 0; 01^ = ''^'t'^ii^) 7^ 0- Evidently the zero-mode solution 0i'(x) is determined up 
to an arbitrary admixture of 0^. When multiplying these linear equations by one 
easily proves with the help of Eqs. (fT^ . (fTS|) that 



[h- - \,)<PUx) = {h- -\,)<Pt{x) = C<Pt{x)- S-=f ^ 1, (20) 




where C is an arbitrary real constant. If Ai 7^ A2 then by the redefinition 
(Ai — A2)0i' = (Ai — A2)0i" + C0^(x) one arrives at the canonical diagonal form S~. 
However in the confluent case, Ai = A2 = A, C 7^ it is impossible to diagonalize and 
by a proper normalization of the zero-mode (pf one gets another canonical form S~ of 
this matrix - the elementary Jordan cell |112j . 

(/i--A)0+(x) = O; {h- -X)<l>t{x) = <j>t{xy, S-=(^^^^. (21) 

We display it to emphasize that in the confiuent case the zero-mode (pf is not anymore 
a solution of the Schrodinger equation but it is a so-called adjoint solution |113j which 
can be obtained by differentiation, (pf = d(j)2 /d\ . Yet the intermediate Hamiltonian 
h = r^r^ + X = Tg^r^ + A is well defined and therefore the intermediate isospectral 
partner 4>t{x) of the zero- mode 4>t{x) is a solution of Schrodinger equation with the 
above Hamiltonian. The analysis of the matrix is similar. Thus we have established 
that in general the Hamiltonian projection onto the subspace of (hermitian) supercharge 
zero-modes is not diagonalizable but can be always transformed into a canonical Jordan 
form. 

To accomplish the description of Polynomial SUSY algebras generated by a 
second-order ladder one should take into consideration also the degenerate case when 
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Ai = A2 = A, C = 0. For this choice the matrix S~ is automatically diagonal and both 
zero modes (f)i2{x) are (independent) solutions of the Schrodinger equation with the 
Hamiltonian h~. Then it can be proved IHS] that the intertwining operator is 
just a linear function of this Hamiltonian, q2 = X — h~. Hence the intertwining is 
trivial h~ = and such supercharges must be eliminated. For higher-order SUSY the 
removal of such blocks in supercharges may lead to ladder irreducible SUSY algebras 
(see Section 4). 

The very fact that the Hamiltonian is represented by finite matrices is 
interpreted sometimes [7S1 ES] as a phenomenon of quasi-exact solvability (QES) 
|114| I115j : this point needs a certain comment. If one is seeking for some formal 
solutions of the Schrodinger equation, not necessarily normalizable and regular then 
such a formal QES can be accepted. But for the spectral design we impose physical 
boundary conditions and requirements of normalizability which are essential to define 
the energy levels properly. Then QES for physical wave functions is achieved only if 
one or both of eigenvalues belong to the energy spectrum of the Super-Hamiltonian. 
Obviously it is an exceptional situation which is not granted by the Polynomial SUSY 
itself. 

Let us complete this section with the general description of the A^-step ladder which 
entails the Polynomial Superalgebra of A^th-order. We introduce a set of first-order 
differential operators for intermediate intertwinings , 

r^ = Td + Xi{x), / = l,...,iV, (22) 

and the relevant number of intermediate super potentials Xi{^)- The set of the initial, 
= ho, the final, h~ = hi\f and intermediate Hamiltonians, hi = —d'^ + vi{x) consists of 
Schrodinger operators, so far nonsingular and real ones. They obey the ladder relations 
("gluing"), 

hi = rf ■ rf + \i = rf_^-^ ■ r^" ^ + Az+i, / = 1, . . . , - 1, 

h^ = h' = ■ rjf + X^, ho = h'^ = ■ + Ai. (23) 

These gluing relations are provided by the (dressing) chain equations on super- 
potentials, 

vi{x) = ixi{x)f + ixiix))' + \i = {xi+i{x)f - ixi+iix))' + (24) 

The corresponding intertwining (Darboux) transformations hold for each adjacent pair 
of Hamiltonians, 

hi-i ■ = rl ■ hi, rf ■ hi^i = hi ■ rf , (25) 
and therefore the chain of overlapping SUSY systems is properly built, 

[Hi, Ri] = [Hi, R]] = 0, Hi-Xi = {Ri, Rj}, (26) 

This Chain Supersymmetry can be equally converted into a A^-order Parasupersymmetry 
similar to Eqs. (fT^ which however we do not need for the further construction. 
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Now let us disregard a chain of intermediate Hamiltonians between and h~ 
and proceed to the Higher-derivative ~ Polynomial ~ Nonlinear SUSY algebra for the 
Super-Hamiltonian H given in Eq.(^. The intertwining between and h~ is realized 
by the Crum-Darboux operators, 

g+=r+...r^, q-=r^...r^. (27) 

The SUSY symmetry [H,Qn] = [iJ, Qjv] = 0, is still performed by the supercharges 
of the same matrix structure (fTB|) and the Super-Hamiltonian is represented by finite- 
dimensional matrices on the subspaces of supercharge zero-modes, 

N 

qUf{x) = 0; t = l,2,...,N; (x) = S^, <Pf (x) , (28) 

i=i 

due to intertwining relations (j21). In terms of the constant matrices , the algebraic 
closure is just represented by a non-linear SUSY relation |H3 [73] , 

N 

^Vn{H) = 1[{H-\i). (29) 
(=1 

The disposition of the real roots A; against the energy levels of the Super-Hamiltonian 
is assumed to provide the positivity of the superalgebra (see next Section). Namely 
if m lowest energy levels, {Ej}, j = 0, . . . , m — 1 are among the roots, {Ej} C {A;} then 
all Xi < Em- Again both matrices S''' have the same set of eigenvalues which for the 
ladder construction (j2Sl) consists of Ai, . . . Aat . If the degenerate roots appear then the 
canonical forms S''' of the latter matrices are not necessarily diagonal and may consist 
of Jordan cells. If all intermediate hi are hermitian, nonsingular and superpotentials are 
taken real, then A/ are real and each ladder step is well defined. What will happen if we 
extend the class of Polynomial SUSY algebras admitting complex A; and singular hi ? 

3. Algebraic classification of Polynomial SUSY QM and its functional 
realization: irreducibility of type I, II, and III 

Let us now examine which circumstances may obstruct the SUSY ladder 
decomposition of a Polynomial SUSY algebra. In fact, all elementary "bricks" 
irreducible to a chain of one-step Darboux intertwinings are well revealed for the second- 
order SUSY algebra described in the previous section. 

For a supercharge of second order in derivatives with real coefficient functions one 
can find real zero- modes of the intertwining operators and, further on, the 2x2 matrix 
representation (fTBj) for the Super-Hamiltonian components by matrices S"*^. The 
latter matrices are real but, in general, not symmetric. Therefore the first obstruction 
for the ladder decomposition may arise because the reduction to a Jordan form has not 
given real eigenvalues. For instance, if h'^(j)~{x) = Xeik(l>j^{x) then the eigenvalues of 
S"*" = Ae are imaginary, mutually conjugated ±iX. The possibility of complex pairs of 
mutually conjugated roots in a Polynomial SUSY algebra can be easily read off from 
its closure (fTIHl as for supercharges with real coefficients polynomials V2{H) possess real 



det 



El-S- 



det 



El-S- 
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coefficients. We call this kind of irreducibility to be of type I. Its elementary block 
corresponds to the polynomial V2^\h) = (H + a)^ + d, d > and its analytical 
properties have been investigated in [Slj, |Hlj- Some examples of related isospectral 
potentials are described in |H2]- 

Next, one has to ensure the positivity of the SUSY algebra relation (fTn|l in a 
particular differential realization of a Super-Hamiltonian H with real non-singular 
potentials and the supercharges Qn,Qn (with = 2 in our case) made of differential 
operators with real coefficients. Let the energy spectrum Ej; j = 0,1, Ej < Ej^i of 
H be discrete, for simplicity. Then, 

Vn{E,) = (Q^^.IQjv^,) + (Q]v^,|Qjv*,) > 0, (30) 

if the action of supercharges is well defined in the Hilbert space spanned by 
eigenfunctions of a Super-Hamiltonian. It can be easily extended on a continuum energy 
spectrum as well using wave packets. 

Thus for regular potentials the allowed disposition of polynomial roots <^==^ zero- 
modes of a supercharge — provides non- negative values of V{E) for each energy level 
of a Hamiltonian. To be definite, one may have the following cases for the allocation of 
polynomial real roots (for a pair of complex, mutually conjugated roots the positivity 
is obvious) . 

Case A. Ai < A2 < Eq or Ai = Eq; A2 = Ei. 

The related SUSY algebra has a chain or ladder realization. In other words, it is 
reducible, in principle, because one gradually can add/remove Ai and then A2 without 
breaking the positivity of intermediate SUSY algebra. The coincidence of roots and 
energies correspond in the spectral design to deleting/inserting energy levels. For 
instance, if Xi = E^; X2 = Ei then two pairs of zero- modes of can be chosen as 
solutions of two Schrodinger equations with Hamiltonians h^. Repeating the arguments 
of Section 1 one can conclude [201 that the energy levels Eq, Ei may well appear in any 
of the Hamiltonians but each level only once, either in or in . Thus the entire 
variety of spectral design tools happens to be at our disposal: namely, one may delete 
two lowest levels, replace the ground-state level by a different one and add two more 
levels below the ground-state one. 

Case B. < Ai < A2 < Ei or Ej < Xi < X2 < -Ej+i, 1 < j- 

A pair of real roots is placed between adjacent energy levels. If one of the solutions 
with eigenvalues Ai^2 (zero-modes of the supercharges) is normalizable we perform 
the insertion/deletion of an excited energy state. Thus with these means one can 
delete two adjacent excited levels, shift the position of an excited level and add two 
more excited levels between two adjacent ones. Evidently, such an algebra cannot be 
safely decomposed into a chain of two linear SUSY as the removal of any of roots Ai_2 
immediately breaks the positivity in (jHUj) . Then the intermediate Hamiltonian acquires 
inevitably a real but singular potential leading to the loss of isospectrality. The related 
Darboux transformations had been known in SOties ^1]. We call this irreducibility to 
be of type II. Examples and certain theorems are given in [HI], |H2]- 
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Case C. Ej < Ai = A2 < ^j+i, < j. 

This is a confluent case which seems to arise as a hmit of the previous one. However, 
let us remind that the one-dimensional QM does not allow degenerate levels. Besides, 
let's assume that the matrix representation for the corresponding Super-Hamiltonian 
contains a non-trivial Jordan cell. Then the limit becomes quite delicate as one of 
the zero-modes is not a solution of the Schrodinger equation but represents an adjoint 
function |113j (see Eq. (j?H) and discussion afterwards). This is why we specify this case 
as a separate one named to be of type III. With such an intertwining operation one may 
insert /delete odd number of excited levels in an economical way. One may find more 
information on the analytical properties of related potentials in [HZj. 

One can apply these second-order blocks and build an A^th-order Polynomial SUSY 
system. Their general form is again given by Eq. (j29|) allowing the presence of complex 
conjugated roots A^. Let us rewrite it taking into account the possibility of complex and 
degenerate roots, 

n m 

1=1 i=i 

n m 

N = Y.ui + 2Y,fi„ dj>0. (31) 

1=1 j=i 

We stress that in the general case the Hamiltonian projections onto the zero-mode spaces 
of intertwining operators are finite N x N matrices S"^ and the Polynomial SUSY 
algebra can be represented by Eqs. (j29|) . Inequality (j30|l is certainly vahd for Vn{H). 

Irreducible elements of type II are not straightforwardly seen in the structure 
of the Polynomial SUSY algebra and can be unraveled only after the inspection of 
disposition of polynomial roots in respect to energy levels. They fill the chain of 
intertwining operators being even order in derivatives and placing a pair (or few pairs) 
of real roots A/ < ■ ■ ■ < A/+2fc-i (supercharge zero-modes) between two successive 
energy levels Ej < \i < ■ ■ ■ < \i+2k-i < Ej^r', r — 1 < 2k if the intermediate levels 
{Ej+i, . . . , Ej+r~i} C {Xi, . . . , Xi+2k-i}', r - 1 < 2k. The eigenvalues Ej,Ej+r are 
assumed not to coincide with any of polynomial roots. Then the polynomial 

2k-l 

t^SUh) = IK^ - ^'+0; -pSU^m) > 0. (32) 

When the related zero-modes coincide with some eigenfunctions of the Super- 
Hamiltonian the pertinent supercharges create or annihilate particular excited states 
in the components of the Super-Hamiltonian. 

Irreducible elements of type HI fill the chain of Darboux transformations being 
represented by even-order intertwining operators responsible for allocation of even 
number of real confluent roots A/ = A;+i = ■ ■ ■ = A;+2!/;-i (= supercharge zero-modes) 
between two adjacent energy levels Ej < A/ < -Ej+i for some < j, 

vS'im = {H ~ X,r- VS^]{EJ>0. (33) 
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No more than two zero-modes may be solutions of the Super-Schrodinger equation, in 
particular, eigenfunctions of the Super-Hamiltonian if Ej = A/. Other zero- modes are 
adjoint functions |llHj to a solution of the Schrodinger equation. 

Finally, in general, the polynomial in Eq. ()31|) can be factorized into the product of 
the elements fjH^ and (jHHjl with roots located between two successive or adjacent levels. 
The related ladder of Darboux transformations consists of reducible steps as well as of 
few irreducible elements of type I, II, III displayed in Eqs. (jHTj) . (jH^ . ((221) • 

Yet the open question remains whether any irreducible element of type II or III ( ()H2j) 
or (|nSl)) can decomposed into the ladder of second-order irreducible blocks with regular 
hermitian intermediate Hamiltonians between them in the ladder. We are informed that 
essential progress in this direction has been made by A. V. Sokolov and hope to see it 
published soon. 

On the other hand, an experienced SUSY designer may be somewhat puzzled with 
the very existence of irreducible super-transformations. Indeed it is quite conceivable 
that a pair of supercharge zero-modes or even a pair of new excited energy levels of the 
Super-Hamiltonian can be inserted by successive application of first-order intertwining 
(super) transformations between regular Hamiltonians following the ladder algorithm 
described in the previous Section. But the order of the relevant ladder of first-order 
transformations and respectively of the final Polynomial SUSY will be evidently higher 
than two. We come to the problem of possible relationship between first-order reducible 
and irreducible SUSY algebras having the same Super-Hamiltonian. 

The related important question concerns the degenerate roots. By general 
arguments these roots are distributed between different Jordan cells in the canonical 
forms of the matrices S^. One can inquire on how many Jordan cells may coexist 
and if several cells appear then what is their role in the supercharge structure. All these 
problems are clarified with the help of the Strip-off theorem |85j . 

4. From reducible SUSY to irreducible one when equipped by the Strip-off 
theorem 

Let first elucidate the possible redundancy in supercharges which can be eliminated 
without any changes in the Super-Hamiltonian [i.e. preserving the same potentials). 
There exists a trivial possibility when the intertwining operators and p^_^ for N > Ni 
are related by a polynomial factor F{x) depending on the Hamiltonian, 

q^ = F{h^)p%^=p%^F{h^). (34) 

Obviously in this case the appearance of the second supercharge does not result in any 
new restrictions on potentials. 

Thus the problem arises of how to separate the nontrivial part of a supercharge 
and avoid numerous SUSY algebras generated by means of "dressing" It can be 

systematically realized with the help of the following 
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"Strip-off" theorem. 

Let's admit the construction given by Eqs. (j28|l and (j29|) . Then 

a) the matrix (or S+) generated by the Hamiltonians (or h~^) on the subspace 
of zero- modes of the operator (or g^), after reduction into the Jordan form S~ (or 
S^), may contain only one or two Jordan cells with equal eigenvalues A^; 

b) assume that there are n pairs (and no more) of Jordan cells with equal eigenvalues 
and with the sizes z// of a smallest cell in the l-th pair; then this condition is necessary 
and sufficient to ensure for the intertwining operator q^^ (or q^) to be represented in 
the factorized form: 

71 

Q^ = P%.m>^i'h^r', (35) 
1=1 

where p%^ are intertwining operators of order A^^i = N — 2 Yl?=i which cannot be 
decomposed further on in the product similar to with F{x) ^ const. 

Remark. The matrices cannot contain more than two Jordan cells with the same 
eigenvalue A because otherwise the operator \ — would have more than two linearly 
independent zero-modes (not necessarily normalizable) . 
The full proof of this theorem has been performed in |85] . 

Let us illustrate the Theorem in the Example: 
the matrix S~ for the intertwining operator with Jordan cells of different size having 
the same eigenvalue. It is generated by the operators, 

= =F<9 + X) = p'^p'^ + A, qt = —p^p~p~^ = p'^{\ — h^). (36) 



Respectively: 
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As a consequence of the "Strip-off" Theorem one finds that the supercharge 
components cannot be factorized in the form ()H4|1 if the polynomial Vn{x) in the SUSY 
algebra closure does not have the degenerate zeroes. Indeed the SUSY algebra 
closure contains the square of polynomial F{x), for instance, 

q^q+ = F{h-)p],y^F{h-) = F\h-)VNAh-), (38) 

where Vni{x) is a polynomial of lower order, A^i < N — 2. Therefore each zero of the 
polynomial F{x) will produce a double zero in the SUSY algebra provided by (jHEl)- 

Thus the absence of degenerate zeroes is sufficient to have supercharges without 
redundancy in the sense of Eq. (jMjl . However it is not necessary because the degenerate 
zeroes may well arise in the confluent ladder construction giving new pairs of isospectral 
potentials [51] . 

Now we proceed to uncover the origin of irreducible, type-II and -III transformations 
based on the strip-off factorization. For clarity let us consider an example of irreducible 
SUSY of type II with supercharges of second order in derivatives (see previous Section). 
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Suppose that it realizes insertion of two new energy levels between the ground and first 
excited states. Then three lowest energy levels Eq < Ei < E2 are of importance to 
study the relevant SUSY systems: the ground state level is degenerate between SUSY 
partners h'^ and h~, i.e. Eq = Eq whereas the two excited levels are present only in 
the spectrum of h~. 

One can use the ladder construction (PHj) - (f?rj) to prepare the same spectral pattern. 
For this purpose, intertwining operators ()27|) of, at least, fourth order in derivatives must 
be employed. Indeed, one can prescribe the ladder steps for qf as follows: start from 
a pair of isospectral Hamiltonians with ground state energies E^; generate the level Eq 
in the Hamiltonian h'^ using the intertwining operators rf , then sequentially create in 
the spectrum of h~ the state with energy E2 < E^ by means of rf, next the energy 
level El < E2 using rf and finally the ground state with energy Eq < Ei exploiting 
rf. These elementary steps are reflected in zero-modes of the intertwining operators 
(it = rfrtr^rf and = r^r^r2ri . Namely the ground state of is a zero- mode of 
rf {i.e. of gj) whereas the eigenstates of h~ corresponding to Eq, Ei, E2 are annihilated 
by the product r2rtrf {i.e. by ^4') according to Eq. (j^H)). In particular, the ground 
state of h~ is a zero- mode of r^. As the ground state energies coincide for the 
Hamiltonian projections on the qf zero-mode space - the matrices are, in general, 
not diagonalizable but have one rank-two Jordan cell each. Thus, for instance. 



(39) 



where a non-zero constant C can be normalized to C = 1. The canonical Jordan form 
S~ in fl39j) is achieved by means of re-factorization qf = rfrfrfrf = rfrfrfrf so that 
the annihilation of ground state for h~ is associated now with rf- Respectively, the 
Polynomial SUSY algebra shows up one degenerate root, 

V,{H) = {H~ Eq)\H - Ei){H - E2). (40) 

The Strip-off theorem tells us that this fourth-order algebra cannot be optimized to a 
lower-order one because there is no replication of roots in different Jordan cells of 
matrices. However one may perform fine-tuning of Darboux transformation parameters 
to provide the constant C = in ()39|1 . This peculiar choice moves the SUSY system into 
the environment of the Strip-off theorem as now two rank-one cells in fj39p contain the 
same eigenvalue Eq . The SUSY algebra is still given by Eq. (j40|l but the intertwining 
operators reveal a redundancy, 

qf = {EQ-h-)qf. (41) 

By construction the left-hand side of this relation is fully factorizable in elementary 
binomials rf with hermitian nonsingular intermediate Hamiltonians. But in the right- 
hand side the operator qf = rfrf does not admit a further factorization with a 
nonsingular intermediate Hamiltonian because after removal of the redundant factor 



/ ^0 








c \ 





E2 














El 





I 








Eo ) 



( Eq 


C 













Eq 
















E2 







I 








El 


) 



Nonlinear Supersymmetry 



14 



{h — Eq) such a factorization is forbidden by the positivity of the SUSY algebra, 

Eq. dnni)- 

One can easily extrapolate the previous argumentation to the case of additional 
degeneracy of excited levels Ei = E2 to analyze the irreducible SUSY of type III. Thus 
we reach the important conclusions that: 

a) the factorization (jTfj) of intertwining operators is not unique and there exist 
options to have more reducible ladders and less reducible ones with a larger number of 
singular intermediate Hamiltonians; 

b) (some of) irreducible algebras of type-II and -III can be identified with special cases 
of fully reducible ladder-type algebras when the Hamiltonian projections have an 
appropriate number of pairs of Jordan cells with coinciding eigenvalues; 

c) there are many (almost) isospectral systems with different pattern of excited states 
which cannot be interrelated with the help of irreducible Darboux transformations of 
type-II or -III but can be built with the help of higher-order reducible SUSY ladder. 

Yet one may substantially gain effectiveness when the spectral design program 
allows to apply the irreducible transformations of type-II or -III in order to embed a 
couple of energy levels between two excited ones. Thus a more rigorous investigation 
of the relationship between the reducible and irreducible intertwinings is welcome. 
Especially important is the proof that any type-II, -III irreducible SUSY can be 
embedded (for the same Super-Hamiltonian) into a reducible ladder SUSY. 

5. More supercharges <^=^ Extended SUSY <^=^ Hidden symmetry 

The possibility of two supercharges for a quantum SUSY system was mentioned in 
[116^ (see the preprint version). Namely the conserved supercharges Q, with complex 
coefficient functions in intertwining components accounted for two SUSY algebras 
for a hermitian Super-Hamiltonian H: one for their "real" parts K, K'^ and another 
one for their "imaginary" parts P, where real and imaginary parts are referred to 
coefficients in the intertwining operators = + ip%^ . 

Let us examine the general possibility to have several supercharges for the same 
Super-Hamiltonian. First we remind that a number of supercharges can be produced 
with the help of multiplication by a polynomial of the Hamiltonian (see Section 4). 
Certainly such supersymmetries are absolutely equivalent for the purposes of spectral 
design and one must get rid of them. As shown in one can always optimize the 
infinite set of possible supercharges so that no more than two nontrivial supercharges 
remain which are used to generate all other ones by "dressing" with polynomials of the 
Hamiltonians. Thus in one-dimensional QM one may have the TV = 1, 2 SUSY only. 

Correspondingly we consider now a general case when the Super-Hamiltonian H 
admits two supersymmetries with supercharges K and P generated by differential 
intertwining operators of order and A^i respectively, 

[H, K] = [H, P] = [H, Kt] = [H, Pt] = 0. (42) 

The second supercharge P is assumed to be a differential operator of lower order Ni < N. 



Nonlinear Supersymmetry 



15 



To close the algebra one has to include all anti-commutators between supercharges, 
i. e. the full algebra based on two supercharges K and P with real intertwining operators. 
Two supercharges generate two Polynomial SUSY, 

{k,K'^}^Vn{H), {p,P^}^VnAH). (43) 

The closure of the extended, Af — 2 SUSY algebra is given by 

{.-.p.}.,^..(*f. ^_;^,). (44) 

Apparently the components of operators TZ, TV are differential operators of A^^+iVi order 
commuting with the Hamiltonians /i^, hence TZ, are symmetry operators for the 
Super-Hamiltonian. However, in general, they are not polynomials of the Hamiltonians 
/i^ and these symmetries impose certain constraints on potentials. 

All four operators Vn{H), Vmi{H), TZ, TZ^ mutually commute. Moreover the 
hermitian matrix describing this J\f — 2 SUSY, 

Z{H) = ( ^^f ^ p^^^^ ) , det [Z{H)] = VnVn, - 7^7^t = 0, (45) 

is degenerate. Therefore it seems that the two supercharges are not independent and by 
their redefinition (unitary rotation) one might reduce the extended SUSY to an ordinary 
Af = 1 one. However such rotations cannot be global and must use non-polynomial, 
pseudo-differential operators as "parameters". Indeed, the operator components of the 
"central charge" matrix 2(H) have different order in derivatives. Thus, globally the 
extended nonhnear SUSY deals with two sets of supercharges but when they act on 
a given eigenfunction of the Super-Hamiltonian H one could, in principle, perform 
the energy-dependent rotation and eliminate a pair of supercharges. Nevertheless this 
reduction can be possible only after the constraints on potentials have been resolved. 

Let us find the formal relation between the symmetry operators TZ, TZ^ and the 
Super-Hamiltonian. These operators can be decomposed into a hermitian and an anti- 
hermitian parts. 

The operator ;B is a differential operator of even order and therefore it may be a 
polynomial of the Super-Hamiltonian H. But if the operator S docs not vanish 
identically it is a differential operator of odd order and cannot be realized by a polynomial 
of ii". 

The first operator plays essential role in the one-parameter non-uniqueness of the 
SUSY algebra. Indeed, one can always redefine the higher-order supercharge as follows, 

X(^) ^K + CP, {k^^^ , K^'^^^ } = (H) , (47) 
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keeping the same order of Polynomial SUSY for arbitrary real parameter (. From 
dHj) one gets, 

2CB{H) = vli\H) - Vn{H) - e-pNAH), (48) 

thereby the hermitian operator S is a polynomial of the Super-Hamiltonian of the order 
Ni, < N — 1. Let's use it to unravel the Super-Hamiltonian content of the operator S, 

£\H) = Vn{H)VmAH) - B\H), (49) 

which follows directly from and (0^1) • As the (nontrivial) operator S{H) is a 
differential operator of odd order Ne it may have only a realization non-polynomial in H 
being a square root of ()49|1 in an operator sense. This operator is certainly non-trivial if 
the sum of orders + A''i of the operators and is odd and therefore N^, = N + Ni. 
The opposite statement was also shown in '851, namely if the symmetry operator is non- 
zero then for any choice of the operators and an optimal set of independent 
supercharges (possibly of lower orders) can be obtained which is characterized by an 
odd sum of their orders. 

The existence of a nontrivial symmetry operator S commuting with the Super- 
Hamiltonian results in common eigenstates which however are not necessarily physical 
wave functions. In general they may be combinations of two solutions of the Shrodinger 
equation with a given energy, the physical and unphysical ones. But if the symmetry 
operator S is hermitian in respect to the scalar product of the Hilbert space spanned 
by the eigenf unctions of the Super-Hamiltonian H then both operators have a common 
set of physical wave functions. This fact imposes quite rigid conditions on potentials. 

In particular, for intertwining operators with sufficiently smooth coefficient 
functions having constant asymptotics for x — > ±oo the symmetry operator S has the 
similar properties and is evidently hermitian. In this case one has non-singular potentials 
with constant asymptotics and therefore a continuum energy spectrum of H with wave 
functions satisfying the scattering conditions. Thus the incoming and outgoing states, 
ipin{x) and ipout{x), at large x are conventionally represented by combinations of plane 
waves which are solutions of the Schrodinger equation for a free particle, 

'ip{x)\^^_oo — ^ exp{ikinx) + R{kin)exp{-ikinx), 

ip{x)\^^+oo — ^ (1 +T{kout))exp{ikoutx), (50) 

where the reflection, R{kin), and transmission, T{kout), coefficients are introduced. 
Since the symmetry is described by a differential operator of odd order which tends 
to an antisymmetric operator with constant coefficients the eigenstates of this operator 
approach asymptotically individual plane waves ~ exp{±ikx) with opposite eigenvalues 
~ -^kf{k'^) and cannot be superimposed. Hence the eigenstate of the Super-Hamiltonian 
with a given value of the operator £ may characterize only the transmission and cannot 
have any reflection, R{kin) = 0. We conclude that the corresponding partner potentials 
Vi^2 inevitably belong to the class of transparent or reflectionless ones jllYj . Such a 
symmetry may have relation to the Lax method in the soliton theory jTH]- 
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As the symmetry operator S is hermitian its eigenvalues are real but, by 
construction, its coefficients are purely imaginary. Since the wave functions of bound 
states of the system H can be always chosen real we conclude that they must be zero- 
modes of the operator S{H), 

S{H)^, = S{E,)^, = 0, VN{Ei)VNAEi) - B\Ei) = 0, (51) 

which represents the algebraic equation on bound state energies of a system possessing 
two supersymmetries. Among solutions of (j3T|) one reveals also a zero-energy state at 
the bottom of continuum spectrum. On the other hand one could find also the solutions 
which are not associated to any bound state. The very appearance of such unphysical 
solutions is accounted for by the trivial possibility to replicate supercharges by their 
multiplication by the polynomials of the Super-Hamiltonian as discussed in Sec. 4. 



6. A simple but useful example of Extended SUSY 

Let us examine the algebraic structure of the simplest non-linear SUSY with two 
supercharges, 

k^ = d^T 2/(x)9 + h{x) T fix); 

p^=^d + x{x), (52) 

induced by the complex supercharge of second order in derivatives [23 I116j . The 
supersymmetries fl42jl generated by K, and P, prescribe that 

where / are real functions and a, d are real constants. The related superalgebra 
closure for K"^ and P, P^ takes the form, 

{K, = {H + af + d, {P, Pt} = H. (54) 

The compatibility of two supersymmetries is achieved by the following constraint x = 2/ 
and by the nonlinear second-order differential equation 



with solutions parameterized by two integration constants. Therefore the existence of 
two SUSY reduces substantially the class of potentials for which they may appear. 
Evidently Eq. (j^^ can be integrated to, 

(/')' = 4/' + 4a/2 + 4Go/ -d^ $4(/), (56) 
where Gq is a real constant. 
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The solutions of this equation are elhptic functions which can be easily found in 
the implicit form, 

[^^^^ ^£= = ±{x-Xo), (57) 

where /o and xq are real constants. 

It can be shown that they are nonsingular if: 

a) $4(/) has three different real roots and the double root P/2 is either the maximal 
one or a minimal one, 

<^>4(/) = 4(/ - ((/ + - (/3^ - e)) , < e < (3'. (58) 

Then there exits a relation between constants a, d, Gq in terms of coefficients /3, e, 

a = 6-^<0, Go = /3(/32-e), d = (^-^ - ^ . (59) 

The constant /o is taken between the double root and a nearest simple root. 

b) $4(/) has two different real double roots which corresponds in (j^ . to 
Go = 0, = e > 0, a = — e/2, d = — The constant /o ranges between 
the roots. 

The corresponding potentials Vi_2 are well known |117j to be reflect ionless, with one 
bound state at the energy (/3^ — e) and with the continuum spectrum starting from 
Respectively the scattering wave function is proportional to exp{ikx) with k = y/E ~ (3"^. 

In the case a) the potentials coincide in their form and differ only by shift in the 
coordinate ( "Darboux Displacement" [H21), 

ch^ (v^(x - a;|) ' 0) ^v^ + 

and in the case b) one of the potentials can be chosen constant (being a limit of infinite 
displacement), 

2 

ch^ (/5(~-Xo)), 

For these potentials one can elaborate extended SUSY algebra. 

The initial algebra is given by the relations (jSH)- It must be completed by the 
mixed anti-commutators 

{K, Pt} = P}t = B{H) - zS{H), (62) 

where the first term is (see the next section) a polynomial of the Super-Hamiltonian and 
the second one is in general not. In our case the first, polynomial symmetry operator 
turns out to be constant, B{H) = Gq when taking into account (jK^ and Meanwhile 
the second operator reads. 



V,=P', V^2 = /^M 1 - TTTT:;?— TV 1 ' (61) 



S{H) 



Id' - {al + ^Y{x))d-jV'{x) 



(63) 
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in the notations H = —9^1 + V(a;). By construction the operator S{H) reahzes a new 
symmetry for the Super-Hamiltonian. Directly from Eq. (jUHj) one derives that, 

£'{H) = h[{H + af + d]-Gl = iH- E,f{H - P^), (64) 

where Ei, = P"^ — e is the energy of a bound state. Thus (some of) the zero modes of S{H) 
characterize either bound states or zero-energy states in the continuum. However there 
exist also the non-normalizable, unphysical zero-modes corresponding to E = Eh, (3"^. 
We remark that in the case b) only the Hamiltonian h~ has a bound state. Hence 
the physical zero modes of S{H) may be either degenerate [case a), broken SUSY] or 
non-degenerate [case b), unbroken SUSY]. 

The square root of (jHH) can be established unambiguously from the analysis of 
transmission coefficients, 



£{H) = {H-Eh)^H-(3\ (65) 

We notice that the symmetry operator (jHSI), is irreducible, i.e. the binomial (H—Eb) 
cannot be "stripped off". Indeed the elimination of this binomial would convert the 
third-order differential operator ()63p into an essentially nonlocal operator. The sign 
of square root in is fixed from the conventional asymptotics of scattering wave 
functions ~ exp{ikx) and the asymptotics Vi^2 — /^^ by comparison of this relation 
with Eq. 

When taking Eq. into account one finds the mixed anti-commutators of the 
extended SUSY algebra ()62|1 in a non-polynomial form. 



{K, Pt} = {K\ Py = Go- i{H - Eh)^H - (66) 

Thus the "central charge" of this extended SUSY is built of the elements (fHl|) and (jHEI) 
and evidently cannot be diagonalized by a unitary rotation with elements polynomial in 
H. Therefore the algebra must be considered to be extended in the class of differential 
operators of finite order. 

Let us now clarify the non- uniqueness of the higher-order supercharge and its role 
in the classification of the Polynomial SUSY. For arbitrary ( in (jTTjl one obtains 

{K^^\ = H^ + (2a + e)H + a^ + d + 2(^0 = (H + a^f + d^, 

ac = a+^C, d^ = d + 2CGo-ae-^C^-^,{-^), (67) 

where $4(/) is defined in Eq. fl^ . 

For the extended SUSY one can discover that the previous classification (Section 
3) of irreducible ladders (Darboux transformations) may fail. Indeed, the sign of c?^, in 
general, depends on the choice of (. For instance, let us consider the case a) when 

= (C + Pf [(C - Pf - 4(/?^ - e)] . (68) 

Evidently if ( lies between the real roots of the last factor in (j68p then d^ is positive and 
otherwise it is negative. But two real roots always exist because jS"^ > e. Thereby the 
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sign of can be made negative as well as positive without any change in the Super- 
Hamiltonian. Hence in the case when the Polynomial SUSY is an extended one, with 
two sets of supercharges, the irreducibility of type I of a Polynomial SUSY algebra does 
not signify any invariant characteristic of potentials. 

7. Non-stationary Schrodinger equations: intertwining and hidden 
symmetry 

Discussions of symmetry properties of the time dependent Schrodinger equation have 
a long history, see for instance |124j and references quoted therein as well as |125j . In 
these discussions the potentials concerned are mainly time independent, see e.g. |118j . 
|12()j . Here our aim is to elucidate that many of the nonlinear SUSY constructions 
illustrated before can be implemented also in the Schrodinger time dependent framework 



Certainly, general first- and higher-order intertwining relations between non- 
stationary one-dimensional Schrodinger operators can be easily introduced. But already 
in the first-order case the intertwining relations imply some hidden symmetry which in 
turn leads to a specific quantum dynamics when the evolution is described by quantum 
orbits and results in the /^-separation of variables |122j . Second-order intertwining 
operators )122j . |12j] and the corresponding non-linear SUSY give rise to the quantum 
motion governed by the spectrum generating algebras. 

Let us start with intertwining relations of the non- stationary Schrodinger operator 



Here dt = d/dt and dx = d/dx denote the partial derivatives with respect to time 
and position : we will denote these derivatives, if applied to some function /, by 



with, in general, complex-valued functions ^0,^1 and ^2-The possibility of a 
complexification of the intertwining (Darboux) (first and also higher order) operator 
was emphasized by |123j . Note also that in contrast to |119j the formalism of |122j 
allow a priori for a first-order operator in dt- 

For the above defined Schrodinger operator ()69|) the intertwining relation reads 



where the functions {i = 0, 1, 2) and Vi^2 are not independent. It can be also 
represented in the SUSY form Eq. © when the stationary Hamiltonians /i^ are extended 
to the Schrodinger operators iS[V^i^2]5 then 



S[V] =tdt + dl-Vix,t) . 



(69) 




(71) 




(72) 
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Inserting the explicit forms of the Schrodinger operators (j69p and the intertwining 
operator (fTOj) into relation (f7T|) it was found |122j that and may depend only 
on time, i.e. = = The assumption that does not vanish identically then 
leads to the consequence that also the potential difference Vi — V2 does depend only on 
time. This is a rather uninteresting case and, therefore, we set .^0 = without a loss of 
generality. Making now the following choice of appropriate variables ^i{t) = e^^^*^ p{t) 
and (^2(2;, = e^^^''^^ p{t)uj'(x,t) with real /5, positive p and complex to functions one finds 

Vi(x,t) = iu'^(x,t)+uj"(x,t) -iuj(x,t) +a(t) - pit) + ipit) / pit) , 

(73) 

V2(x, t) = uj''^{x,t) - uj"{x,t) - iuj{x,t) + a{t) , 

where a is some time-dependent complex- valued integration constant. Again one may 
set |122j j3 = without loss of generality. Furthermore, one may also take a = because 
it can always be absorbed in u by the shift u ^ u — i J dt a. Hence, we are left with 

V2{x,t) = uj''^{x,t) — uj"{x,t) — iuj{x,t) . 

Here the so-called super-potential u is still not arbitrary as the potentials are to 
be real. This can, for example, be achieved by assuming a stationary real super- 
potential. However it leads to the standard stationary SUSY QM discussed previously. 
Alternatively, we will consider a complex super-potential 

ui^x, t) = ujr{x, t) + iuji{x, t) (75) 
with real functions uji and ui. The reality condition Im Vi = Im V2 = implies 

2(^,)" + p/p = , 2{iUR)\u;iy - M" -uJr = 0, (76) 
which can easily be integrated to 

uji{x,t) = -l^x^ + l-p{t)fi{t)x + 'y{t) , 
^ 4 p{t) 2 

0OR{x,t) = -\np{t) + K[x/p{t)+p{t)) , 

where p and 7 are arbitrary real functions of time and K is an arbitrary real function of 
the variable y = x/ p + p. In terms of these functions the final form of the two partner 
potentials is 

V,,2{x,t) = 4F7,\K'Hy)±K"{y) 



P\t) 



■||x^+(pWMt) + «).-^ + 7(t) (7. 



and the intertwining operator reads 



g+(x, t) = p{t)d, + K'[x/p{t) + p(t)) - [p{t)x - p\t)p{t)) . (79) 
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Let us demonstrate |122j that the non- stationary Schrodinger equation iS[Vi_2]V'i,2 = 
with potentials given in Eq. ()78p (which is equivalent to the intertwining (f7T|) ) admits 
a separation of variables. In fact, after the transformation 

y = x/p{t) + ^^{t) , V^i,2(x, t) = e-^"^("'*Vi,2(z/, t) = ^](x, t)0i,2(z/, t) (80) 

{pit) 

this Schrodinger equation becomes quasi-stationary |126j 

ip'm<Pi,2{y.t) = [-dl + K'\y) ± K"{y)\ <P,^^{y,t) , (81) 

which is obviously separable in y and t. Hence, the solutions of the original Schrodinger 
equations have the general form il){x,t) = Q{y,t)Y{y)T{t) which is known as the R- 
separation of variables |127j . In other words, for any pair of Schrodinger operators 
'5[Vi,2], which admits a first-order intertwining relation (f7T|l there exists a transformation 
(jHUj) to some new coordinate in which the potentials become stationary (see also |121j ). 
Notice that the transformation associated with the special case p{t) = 1 and fi{t) = vt 
with constant velocity v corresponds to the Galileo transformation. See, for example, 
the textbook [T28] . 

This /^-separation of variables is certainly related to the existence of a symmetry 
operator. First, one can directly verify the adjoint intertwining relation for real 
potentials, 

qrS[V,] = S[V2]qr (82) 

where 

q- ^ (g+)t = -p{t)d, + K'(x/p{t) + pit)) + [p{t)x - p\t)p{t)) . (83) 
Then from (f7T|l . (ff^ and (jH^ we obtain the closure of the SUSY algebra, 

[Qt,Q\} = nt, [5^,7^^]=o, g^ = (g^)^ (84) 

where the symmetry operator TZt has the following components 

Rf = qfqf = -p\t)dl + (pp(t){x, d,} - 2p\t)p{t)d, 



2 



K'(x/p{t) + p{t)y\ ± K"(x/p{t) + pit)) + - [p{t)x - p\t)p{t) 



4 

= exp{-^u;,(x, t)} [-dl + K'\y) ± K"{y)\ exp{2a;/(x, t)}. (85) 

Thus the quasi-stationary Hamiltonians in Eq. (jHT|) are just unitary equivalent to the 
symmetry operators Rf. It means that the supersymmetry entails the separation of 
variables because it provides a new symmetry. As a consequence the quantum dynamics 
splits in orbits with a given eigenvalue of the symmetry operator. 



8. Second-order intertwining for stationary potentials: symmetry operators 
and spectrum generating algebra 

Now we will report on the intertwining of a pair of Schrodinger operators iS[Vi] and 
S\y2\ by second-order (intertwining) operators of the form: 

qt{x, t) = G{x, t)dl - 2F{x, t)d^ + B{x, t) . (86) 
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We will explore the connection of the time- dependent SUSY charges with appearance of 
the spectrum generating (oscillator like) algebras for the corresponding Hamiltonians. 

As in the first-order case it can be shown |122j that the inclusion of an additional 
term being of first order in dt leads to the trivial situations where the difference V1 — V2 
depends on the time t only. Furthermore, from the intertwining relation (jTlj) with 
above one can conclude that the function G may not depend on x and similarly to 
the discussion in the previous section it is even possible to exclude a phase. In other 
words, without loss of generality G{x, t) = g(t) and consider from now on an intertwining 
operator of the form 

t) = g{t)dl - 2F{x, t)d^ + B{x, t) . (87) 

In |122j particular solutions of the intertwining relation (ffT|) were constructed with 
g^" as given above . In this section we shall analyze the solutions of the intertwining 
relation (ffT|l for the case where both potentials Vi and V2 are stationary, i.e. do not 
depend on t. 

One class of such solutions is known from [HI]. Assuming a supercharge qf' with 
real coefficient functions independent on t, one finds that the corresponding solutions of 
fl71|) coincide with those of the stationary intertwining relations {—d^ + Vi{x))q~^{x) = 
q+{x){-dl + V2ix)) from P!- 

Here we are interested in more general solutions of !^71\i when operators qt' depend 

on t, 

(idt - h+)q+{x,t) = qt{x,t){idt - h') , (88) 

with standard stationary Hamiltonians = — 2(3;) but explicitly time-dependent 
intertwining operators. 

Let us employ the suitable ansatz with simple t-dependence in (|H7)). 

g+(x, t) = M+{x) + A{t)a+{x) , (89) 

where 

M+{x) = dl-2f{x)d^ + b{x) , a+{x) = + W{x) . (90) 

Here all functions besides A are considered to be real. We also assume A ^ 0. With 
this ansatz the intertwining relation ()88|1 can be shown |122j to yield 

iA = 2m + 2mA , 

h+M+ - M+h- = 2ma+ , (91) 
h^a^ — a^h~ = 2ma^ , 

with real constants m and m. 

We find it interesting to focus on the case m 7^ to explore certain spectrum 
generating algebras. The first equation in (|nT| immediately leads to 

A{t) = moe-^'""' - rh/m (92) 
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with a real mg, and 

g+(x,t) =dl-( 2f(x) + -)d, + b(x) - - W(x) + moe-^'""' a+ (x) . (93) 
V mj m 

It is obvious that without loss of generality we may set m = as a non-vanishing m 
may always be absorbed via a proper redefinition of / and 6, i. e. of the operator M. 

As a consequence, the second relation in ()9H1 leads to a second-order intertwining 
between and . This has already been considered in [S^ and it was found that the 
potentials Vi, V2 and the function h can be expressed in terms of / and two arbitrary 
real constants a and d: 

fix) r\x) d 



Vi,2(x) = t2/'(x) + /2(x) + 



2/(x) 4/2(x) 4/2(x) 



The corresponding second-order SUSY algebra generated by the supercharge M is 
similar to (jMj) . 

{M,M^} = {H + af + d = V2{H). (95) 

Note that cases with < are reducible ones. 

One may find some similarities between the present intertwining algebra (j9T| and 
the extended SUSY relations discussed in Section 5. But we emphasize that now for 
m 7^ the last relation in (^1]) does not generate a second SUSY. Rather it creates the 
equivalence of relatively shifted spectra of two Hamiltonians and h~ which is typical 
for spectrum generating algebras. Specifically 

a^a" = h'^ ~ m + c; a~a^ = h~ + m + c, (96) 

where c is a real constant. Therefore the refiectionless potentials found in Section 5 
are produced only in the limit of m = 0. For this reason we use here the notations for 
relevant operators different from those ones in Section 5. 

The genuine spectrum generating algebra for stationary Hamiltonians /i^ can be 
derived from Eq. (jnTjl 

[h+, G+] = -2mG+ , G+ = M+a- , 
\h+, gU = 2mG\. , G\ = a+M- , 

(97) 

[h-,Gl] = 2mGl , GI = M-a+ , 
[h-,G^] = -2mG_ , G_ = a-M+ , 

where a~ = {a^Y and M~ = {M~^)K The closure of this spectrum generating algebra 
is a polynomial deformation of Heisenberg algebra jl29j . 

[GlG±] = F^{h^). (98) 

The explicit form of the polynomials F^{x) can be obtained with the help of Eqs. (j^Hl 
and for m = 0. For instance, 

GIg+ = {h+ -m + c) V2{h^ - 2m); 

G+Gt = (/i+ + m + c)P2(/i+), (99) 
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where the notations from Eqs. (|95|) and (j96p are employed. The polynomials turn 
out to be different for the isospectral partners h^, 

F^{h-^) = -6m{h^y + 4m(2m - 2a - c)/i+ - 2m[a2 + d + 2{a - m){c - m)]; 

F-{h-) = -Qm{h-f - 4m(2m + 2a + c)/i^ - 2m[a^ + + 2(a + m)(c + m)]. (100) 

Hence the two spectrum generating algebras are, in general, different that is essentially 
due to the shift in intertwining relations (jHSl)- There is a formal discrete symmetry 
between their constants and Hamiltonians , a, c =^ —h~, —a, —c. 

The intertwining relation (jHHjl and its adjoint give rise to the symmetry operators 
Qt~1i~ qiqt for {idt — h~^) and {idt — h~), respectively. Using Eqs. and 
after elimination of polynomials of the Hamiltonians these symmetry operators may 
be reduced to the form 



R+{x,t) = mo [e2^'^*G'+ + e-^^™*^! 



R {x, t) = rriQ 



(101) 



As our potentials do not depend on time the operators R^{x,t + A) with a time shift 
A are also symmetry operators for the same Schrodinger equation, 

'Rtt+AvSt]=0. (102) 

In particular, time derivatives R^{x,t) of hermitian symmetry operators R^{x,t) form 
an independent set of hermitian symmetry operators which do not commute between 
themselves. Similar results have also been obtained in |118j using a different approach. 
We see that the non-stationary SUSY delivers the symmetry operators which encode 
the entire set of spectrum generating algebras 

^2imtQ ^ -^R±(x, t) - —^R^ix, t). (103) 

2mo ^ ' Ammo 
The natural question concerns the reducibility of the second-order intertwining 
operator (jHTjl to a pair of consecutive first-order operators. Progress in classification of 
possible irreducible transformations has been made in |122j though more work must be 
done toward the full classification. 



9. Conclusions and perspectives 

The purpose of this short review has been two-fold: to elucidate the recent progress 
in Nonlinear SUSY realization for a broad community of spectral designers and to 
draw reader's attention to a variety of SUSY extensions which yield new QES potential 
systems and illuminate some old ones. With the experience from the previous sections 
the general SUSY QM can be thought of as governed by the extended nonlinear SUSY 
algebra with J\f pairs of nilpotent supercharges Qj, and a number of hermitian 
hidden-symmetry differential operators Ra = R^, [Ra,R^] =0; < < M. 
Such a SUSY algebra takes the modified form, 

[Ra, Qk] = [Ra, Ql] = 0; {Q„ Qk} = {Q], Ql} = 0; 

{Q„Qi}=V{R^); (104) 
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We notice that, first, the Super-Hamiltonian itself is included into the set of symmetry 
operators, say for a = 0, Rq = H and, second, not all the symmetry operators are 
necessarily present in the algebraic closure ()1()4|1 (see Sections 7,8). 

On the other hand, there remains a plenty of open questions and challenges to be 
solved. 

• In the first half of this paper we have given a systematic analysis of reducibility 
vs. irreducibility of type I, II, III in the one-dimensional QM. However the higher- 
order irreducibility needs more efforts to prove the exhaustive completeness of the 
classification in Section 3. 

• For SUSY extensions it may be of interest to find pairs of (quasi)isospectral 
potentials admitting hidden symmetries which are related by a type-II irreducible 
Darboux transformation. 

• The irreducibility classification for non- stationary potentials as well as the existence 
of extended SUSY is very welcome to be investigated and new applications to be 
found, in particular, to explore spectrum generating algebras (see Section 8). 

• The similarity of the Schrodinger equation to the Fokker-Planck one allows 
jnni I122[ I13()j to find the SUSY scheme for generating new solutions of the latter 
equation. One can be tempted to develop a more exhaustive analysis of how to 
produce SUSY clones using the ideas of the conventional nonlinear SUSY outlined 
here. But attention should be focused on the non-hermiticity of the Fokker-Planck 
operator and on the fact that the equivalent of the wave function is a positive and 
properly normalized probability function. 

• Matrix (coupled channel) systems represent a rich and not fully scanned field of 
extended SUSY systems with hidden symmetries. While certain interesting matrix 
potentials have been explored [IS], |131j - |134j it is clear that in this case the way 
to a comprehensive understanding of irreducible building blocks for spectral design 
is still long. 

• The polynomial SUSY in two dimensions has already brought a number of examples 
of new type of irreducible SUSY with hidden symmetries of higher-order in 
derivatives jSSl EH I135j . One may expect a variety of new types of irreducible 
SUSY for third-order (and higher-order) supercharges as well as new discoveries in 
three dimensions. 

• Complex potentials ( |136j - P^^ ) seem to offer less problematic generalizations of 
QM with hermitian Hamiltonians as compared to matrix and multi-dimensional 
QM. Therefore many of tools and results of one-dimensional SUSY QM are 
expected to be applicable when a potential is complex [H3]. However, as it was 
recently remarked in |141j there exist non-hermitian Hamiltonians which are not 
diagonalizable but at best can be reduced to a Jordan form. For the latter ones 
special care must be taken to derive the isospectrality and to build SUSY ladders. 
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